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ABSTRACT 

We present and discuss the mean rest-frame ultraviolet spectrum for a sample of 81 Lyman Break 
Galaxies (LBGs) selected to be B-band dropouts with a mean redshift of z = 3.9 and apparent 
magnitudes zab < 26. Most of the individual spectra are drawn from our ongoing survey in the 
GOODS fields with the Keck DEIMOS spectrograph described in earlier papers in the series, and we 
have augmented our sample with published data taken with F0RS2 on the VLT. In general we find 
similar trends in the spec tral diagnostic s to th ose found in the earlier, more extensive survey of LBGs 
at z = 3 undertaken by iShaplev et al.l ()2003[ ). Specifically, we find low- ionization absorption lines 
which trace the presence of neutral outflowing gas are weaker in galaxies with stronger Lya emission, 
bluer UV spectral slopes, lower stellar masses, lower UV luminosities, and smaller half-light radii. 
This is consistent with a physical picture whereby star formation drives outflows of neutral gas which 
scatters Lya and gives rise to strong low-ionization absorption lines, while increasing the stellar mass, 
size, metallicity, and dust content of galaxies. Typical galaxies are thus expected to have stronger 
Lya emission and weaker low-ionization absorption at earlier times (higher redshifts) . Indeed, our 
mean spectrum at z ~ 4 shows somewhat weaker low-ionization absorption lines than at z — S and 
available data at high redshift demonstrate that this evolutionary trend continues. Although the 
total absorption by low-ionization transitions weakens at high redshift, the fine structure emission 
lines are stronger suggesting a greater concentration of neutral gas at small galactocentric radius 
(< 5 kpc). In conjunction with earlier results from our spectroscopic survey which demonstrated 
an increased fraction of LBGs with Lya emission at higher redshift, we argue that the reduced low- 
ionization absorption is likely caused by a decrease in the covering fraction and/or velocity range 
of outfiowing neutral gas at earlier epochs. At present we cannot distinguish between differences in 
the covering fraction, outfiow kinematics and geometry as the underlying cause of this interesting 
evolution. However, our continuing survey will enable us to extend these diagnostics more reliably 
to higher redshift and determine the implications for the escape fraction of ionizing photons which 
governs the role of early galaxies in cosmic reionization. 
Subject headings: galaxies: high redshift - galaxies: evolution 



1. INTRODUCTION 

Considerable progress has been made over the past 
decade in charting the demographics of high redshift 
galaxies. Multi- wavelength surveys have defined the lu- 
minosity functio ns of UV and sub-mm selected star- 
forming sources (IReddv fc Steidell[2009t IWardlow et all 
1201 It) as well as th e coeval population of quiescent mas- 
sive red galaxies (jBrammer et al.|[201lD . Spitzer data 
has revealed the time-dependent stellar mass density - a 
complementary quantity which repres ents the integral o f 
the past star formation activity (e.g. iStark et al.l 120091 ). 
Through these surveys a well-defined picture of the his- 
tory of star formation and mass ass embly over < z < 
6 has been empir ically determin ed ("Hop kins fc BeacomI 
[2001 Ellis 2001 IRobertson eta l. 2010). The redshift 
range 2 < z < 3 corresponds to the peak of star forma- 
tion activity where the Hubble sequence starts to emerge, 
and the earlier era corresponding to 3.5 < z < 5 is an 
even more formative one where mass assembly was par- 
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ticularly rapid. 

Intermediate dispersion spectroscopy of carefully- 
selected Lyman break galaxies (LBGs) has been particu- 
larly important in defining population trends that cannot 
be identified from photometric data alone . A very influ- 
ential study at z ~ 3 was undertaken by IShaplev et al.l 
([2003) who used composite Keck LRIS spectra of vari- 
ous subsets of nearly 1000 LBGs to examine the role of 
hot stars, Hii regions and dust obscuration, as well as to 
measure the outfiow kinematics and absorption line prop- 
erties of neutral and ionized gas. Composite spectra are 
particularly useful for measuring weak lines which cannot 
be studied in detail for individual objects. Through these 
careful studies, a detailed picture of the mass-depen dent 
evolution of LBGs has emerged (see IShaplev! 1201 ll for a 
recent review). 

In earlier pa pers in this series (|Stark et al.ll20ia here- 
after Paper I; IStark et al.l 1201 ll hereafter Paper II), we 
introduced an equivalent spectroscopic survey of LBGs 
selected from a photometric catalog of more distant 
LBGs with 3 < z < 7 in the Gre at Observatories 
Origins Deep Survey (G OODS) fields (iGiavalisco et al. 
'2004!; IStark et al.l 120091) . Whereas the IShaplev et al. 
(2003) study targetted the study of LBGs close to the 
peak of activity in the overall cosmic star formation his- 
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tory, this earlier period corresponds to a less well-studied 
era when the rate of mass assem bly is particu l arly r apid. 
From photometric data alone, iStark et all (|2009[ ) de- 
duced some significant changes in the characteristics of 
star formation at z ~ 4 — 6 compared to later times, for 
example a shorter timescale of activity (~ 300 Myr). We 
considered it crucial to understand these changes in LBG 
properties if these galaxies are to be used as probes of 
cosmic reionization at higher redshifts. 

At the time of writing, our Keck survey is continuing 
with increasing emphasis at high redshift. Paper I pre- 
sented the first substantial results from a survey of LBGs 
at 3.5 < z < 6 observed with the Keck/DEIMOS spec- 
trograph. Paper II augmented this data with a further 
sample following more ambitious exposures focusing pri- 
marily on 2 ~ 6 LBGs. Incorporating a sample of ESO 
VLT spectra, retrospectively selected using similar pho- 
tometric criteri a as those for the Keck san i ple from th e 
F0RS2 study of I Vanzella eFaH (l2005l [20061 [20081 12001) . 
the current dataset amounts to a sample of 546 galaxies 
over the redshift range 3.5 < z < 6.3. 

Our earlier papers in this series concentrated primarily 
on the rate of occurrence of Lyman a (Lya ) emission in 
our spectra (the "Lya fraction" ) . The overall goal was to 
understand the significantly different evolutionary trends 
in the luminosity functions of LBGs and narr ow-band 
selected Lya emitters (LAEs. [Ouchi et aI1l200 8^ prior to 
the use of the Lya fraction a s a test of when reionization 
ended (jSchenker et al."2011') . Paper I confirmed a result 
found bv .Shaplev et a l. (2003) at lower redshift, namely 
that Lyo; emission is more frequent in lower luminosity 
LBGs rising to a high proportion ~ 50% at Muv = —19. 
More importantly, the Lya fraction was found to rise 
modestly with redshift over 3 < z < 6. By correlating 
the visibility of Lya emission with UV continuum slopes 
derived from the HST photometry, it was argued that 
these trends in the visibility of line emission most prob- 
ably arise from different amounts of dust obscuration. 
Reduced dust extinction in lower luminosity LBGs and 
those at higher redshift has also be en deduced from stud- 
ies of larger photo metric samp les (H ouwens eraI1[2009[ 
see also Reddv fc Steidelll200"l . 

Paper I also discussed the possibility that the cover- 
ing fraction of hydrogen may be lower in low-luminosity 
LBGs. Strong Lya emission in luminous LBGs is often 
associated with low equivalent width interstellar absorp- 
tion lines arising fro m a non-uniform covering fraction 
of ne utral hydrogen ([Quider et al.l 120091 : IShaplev et al.l 
l2003f ). This trend suggests that the high Lya fraction 
in faint LBGs is partially due to a lower covering frac- 
tion, which would imply that Lyman continuum photons 
may more easily escape from intrinsically faint z ~ 4 — 6 
galaxies. Such a result would have great importance in 
understanding the r ole oi z > 7 galaxies in maintaining 
cosmic reionization ([Robertson et al.ll20100 . 

The present paper represents our first analysis of the 
spectral properties of z ~ 4 — 5 LBGs derived from com- 
posite spectra in th e manner pioneered at z ~ 3 by 
IShaplev et al.l (|2003f ). The large database now amassed 
following the campaigns at Keck and the VLT makes a 
similar study now practical in the redshift range where 
there is evidence of increased short-term star formation 
and the mass assembly rate is particularly rapid. Via de- 
tailed studies of low-ionization absorption line and emis- 



sion line profiles, we aim to examine possible changes 
in the kinematics and covering fraction of neutral gas, 
which affects the strength of Lya and the escape fraction 
of ionizing photons, as we approach the reionization era. 
As our redshift survey continues, in this paper we focus 
on a sample of galaxies with 3.5 < z < 4.5, selected as 
B-dropout LBGs. Combining our Keck sample with data 
from the VLT (see Paper I for details), herein we exam- 
ine the spectral features and trends in composite spectra 
drawn from a sample of 131 galaxies. 

A plan of the paper follows. We briefiy review the 
spectroscopic observations and their data reduction in 
§2; much of the relevant discussion is contained in Paper 
I. In §3 we describe the selection of individual spectra 
that we consider appropriate for forming the composite 
mean spectrum at z ~ 4 and the associated sources of un- 
certainty. §4 examines the mean spectrum in detail and 
introduces the various diagnostic features in the context 
of a physical model for LBGs of different masses and star- 
formation rates. In §5 we compare spectroscopic trends 
grouped by observable properties such as mass and lu- 
minosity with those found at z ~ 3 by IShaplev et al.l 
(2003). In §6 we discuss those trends which appear to be 
redshift-dependent, discussing implications for the role 
of early star-forming galaxies in cosmic reionization. Fi- 
nally we summarize our results in §7. 

Throughout this paper, we adopt a fiat ACDM cos- 
mology with f^A — 0.7, Hm — 0.3, and Hq = 70 /lyo 
kms~^Mpc~^. All magnitudes in this paper are quoted 
in the AB system fOke.1974 .). 

2. OBSERVATIONS AND DATA REDUCTION 

The rationale and procedures used to undertake our 
spectroscopic survey of LBGs over 3 < z < 7 were in- 
troduced in detail in Paper I and the interested reader 
is referred to that paper for further detail. Here we 
recount only the basic details. Our target LBGs were 
selected as ^ or i'-band "dropouts" based on deep 
photometry in the two GOODS fields. The photometric 
catalog used in our analysis will be described in detail 
in Stark et al. (2011, in preparation). The selection 
and photo metric appr o ach is largely similar to that de- 
scribed in IStark et "all ([2009D and Paper I, but we high- 
light two key updates. First the selection is performed 
on the v2 GOODS ACS catalogs. Second, we utilize 
deep ground-based near-IR imag ing in GOODS-N ob- 
tained from WIRCAM on CFHT ([Wang et al|[2010l ) and 
deep HST Wide Fie ld Camera 3 / IR imaging of GOODS- 
S from CANDELS (IGrogin et alllMTl [Koekemoer et all 

[mil). 

2.1. Keck/DEIMOS 

The majority of spectra we present are taken from 
an ongoing survey with t he DEep Ima ging Multi-Object 
Spectrograph (DEIMOS: [Faber et al.li2003.' ) on the Keck 
II telescope. In this paper we have used observations 
taken in 2008 April and 2009 March (masks GN081, 
GN082, GN083, GN094, and GN095 in Paper I), which 
targeted a total of 261 B-drop and 88 t^-drop galaxies. 
These data were taken with the 600 lines mni^^ grating, 
covering the wavelength range 4850-10150 A with a res- 
olution of ~ 3.5 A. We do not use data at wavelengths 
A > 9200 A which are affected by strong and variable 
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absorption by atmospheric water vapor. The seeing was 
typically O'.'S and ranged between 0'.'5 and I'.'O during the 
observations. 

All data were reduced and calibrated using a mod- 
ified version of the IDL pipeline SPEC2d, developed 
specifically for D EIMOS by the DEEP2 survey team 
(|Davis et al.ll2003D . The data were reduced as described 
in Paper I, with the addition of two important modifi- 
cations. First, the continuum traces of all target galax- 
ies and other objects occupying the same slit were care- 
fully masked to exclude object flux from the sky back- 
ground model. Second, the b-spline fit to the sky back- 
ground was modified to include a 2nd-order polynomial 
fit to the spatial dimension. These modifications sig- 
nificantly improved the sky subtraction, particularly at 
long wavengths A > 7000 A where bright sky lines can 
be problematic. 

The reduced one- and two-dimensional spectra were 
visually inspected using the IDL program SpecPro 
([Masters &: Capa"^ l2011[ ). Spectra which suffered from 
poor data quality were excluded from further analysis. 
The exclusions included unacceptable amounts of scat- 
tered light within the detector mosaic, defective CCD 
columns, contamination from bright nearby sources, or 
poor sky subtraction, for example arising from the lo- 
cation of the slit on the detector. A few low-redshift 
interlopers were also identified and excluded from fur- 
ther analysis. For the remaining spectra, redshifts were 
measured from either Lya emission (where present) or in- 
terstellar absorption lines. Galaxies identified as hosting 
strong AGN based on the presence of CiV or other strong 
emission lines were excluded. We identified 5 stars, 2 
dusty low-redshift galaxies {z ~ 0.5), 4 AGN, and 134 
star- forming galaxies with secure redshifts z > 3. Af- 
ter rejecting poor data, the final sample consists of 94 
high-quality spectra with accurate redshifts. Examples 
of high redshift dropout spectra are shown in Figure [T] 

2.2. Archival VLT/F0RS2 Spectroscopy in GOODS-S 

To augment our sample of high redshift spectra, we 
have also made use o f data from the FQ RS2 program of 
IVanzella etJaD (|2005l 120061 12008L |2009( ) which targeted 
dropouts in the GOODS-S field. The characteristics of 
that survey in terms of resolution and spectral cover- 
age are very similar to that undertaken at Keck with 
DEIMOS and details can be found in Paper I. Using the 
coordinates provided in the published F0RS2 database 
we queried the version 2.0 ACS catalogs for GOODS- 
S and undertook our own photometric measures and 
dropout selection criteria in an identical fashion to that 
used for our Keck survey. The magnitude distribution of 
the F0RS2 sample is generally weighted towards sources 
brighter than those in the overall Keck survey, but for 
the purposes of constructing the mean spectra discussed 
in this paper, the bulk of the individual spectra are of 
comparable brightness. 

2.3. Redshift measurements 

Care is needed in deriving accurate systemic redshifts 
from rest-frame UV spectra since the strongest features 
trace the kinematics of outflowing gas rather than that 
of the stars. Stellar absorption lines are usually too faint 
to be measured precisely given the signal to noise of the 



spectra. Typically the only features detected in indi- 
vidual spectra are Lya and strong interstellar absorption 
lines such as SiiiA1260, OiA1302-hSiiiA1304, CiiA1334, 
SiivAA1393,1402, SiiiA1526, and CivAA 1548,1550 (Fig- 
ure [T|). Absorption by outflowing gas results in a 
blueshift of interstellar absorption lines. Outflowing neu- 
tral hydrogen along the line of sight leads to a Lya proflle 
which displays very broad (> 1000 kms~^) blueshifted 
absorption and net redshifted emission. The high- 
ionization Siiv and Civ lines also arise in P-Cygni stellar 
winds with broad blueshifted absorption. The magnitude 
of these offsets has been well-quantified for star-forming 
galaxies at z ~ 2.3, and is typically —200 kms~^for inter- 
stellar absorption li nes and -1-500 kms^^for Lya emission 
(ISteidel et al.llMol) . 

To determine accurate redshifts for making compos- 
ite spectra, we restrict our sample to those with red- 
shifts measured from either Lya emission {z^ya ) or low- 
ionization interstellar absorption lines (z/s). Although, 
as discussed above, these are not at the systemic red- 
shift, we follow well-e stablished techniqu es to correct 
for the typical offsets C Steidel et 311120101 ). We do not 
use redshifts based on Lya absorption or high- ionization 
lines (Siiv and Civ) because of the complex and variable 
blueshifts of these features with respect to the systemic 
redshift. We define ZLya as the centroid of the emission 
line and consider only spectra in which the line is de- 
tected at > 5(7 significance. Interstellar absorption line 
redshifts require careful treatment to avoid spurious iden- 
tification of sky line residuals. We consider only the low- 
ionization SiiiA1260, OiA1302-hSiiiA1304, and CiiA1334 
features which are typically found in the highest signal- 
to-noise regions of our spectra. Absorption features at 
longer wavelengths are less reliable due to the higher den- 
sity of strong night sky lines, while shorter wavelength 
transitions are lost in the Lya forest. To measure ab- 
sorption line redshifts, we first estimate the redshift from 
Lya (either in emission or absorption) . We then fit Gaus- 
sian profiles to the spectrum near the expected position 
of the three features above. We require the best fit of all 
three lines to be consistent to within ±500 km s^^ with 
a combined significance of > 5ct. If these conditions are 
met, we define zjs as the weighted mean redshift of the 
three interstellar features. 

With the criteria above, suitable redshifts are available 
for a total of 131 high-quality DEIMOS and F0RS2 spec- 
tra. 91 redshifts are based on measures of ZLya emission 
only, 31 have zjs only, and 9 have both ZLya and z/s- 
The distribution of redshifts and absolute UV magni- 
tudes for this sample is shown in Figure [2l Unless stated 
otherwise, further analysis in this paper is restricted to 
the 81 sources with redshift z < 4.5 and apparent magni- 
tude z'^^ < 26.0. We applied this additional magnitude 
criterion in order to ensure a well-defined continuum sig- 
nal/noise in each individual spectrum. 53 of these 81 
sources are drawn from the Keck survey and 28 from 
F0RS2. 

2.4. Sample bias 

While the sample is 90% complete for dropout selected 
galaxies to an apparent magnitude zab = 25, at fainter 
magnitudes Lya emission or strong interstellar absorp- 
tion features are required for reliable redshifts. This re- 
sults in a bias towards stronger low-ionization absorp- 
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FIG. 1: Examples of individual DEIMOS spectra for 3.5 < z < 4.5 LBGs. From top to bottom: a galaxy at 2 = 3.73 with 
Lya emission and weak interstellar absorption features; a galaxy at z = 3.906 with Lya and interstellar absorption; a faint galaxy 
at 2 = 4.07 with Lya emission and no detectable absorption lines; and galaxy ai z — 3.778 with Lya emission whose spectrum 
is contaminated by scattered light and detector defects. The top three spectra are included in stacking analyses while the 
bottom spectrum is excluded based on poor data quality. The continuum signal-to-noise in each spectrum degrades noticeably 
at wavelengths A > 7000 A due to increased OH sky emission. 



tion lines for the fraction (~ 35%) without Lya . We can 
quantify this bias through the detectability of the aver- 
age low-ionization absorption line strength shown in the 
composite spectrum in Figure [3] (see fjS]). We find that 
we can measure redshifts for this average absorption line 
strength with 90% completeness at zab — 24.7, declining 
through 50% at zab = 25.2 to zero at zab > 25.5. Of 
the galaxies with detected absorption line redshifts, those 
with ZAB > 25.2 have line strengths only 10% stronger 
than for those with zab < 25.2. Clearly this is a small 
effect. 

The bias toward stronger Lya emission for fainter 
galaxies is manifest in Figure[2]where a paucity of objects 
with WLya < 20 A can be seen at faint magnitudes. This 
bias was fully quantified in Paper I using Monte Carlo 
simulations. For galaxies at 3.5 < z < 4.5, a sample 
completeness of 95% is reached at WLya ?i 20 A for 
z'^g = 26, and WLya > 7 A for z'^g = 25. For the 
sample presented in this paper, the least biased subset 
is that with strong Lya emission, followed by that with 
bright apparent magnitudes z'^^^ < 25. 

A final issue in considering composite spectra is that 



these are comprised of individual spectra across 3.5 < 
z < 4.5 with different rest-frame wavelength ranges. 
As our spectra generally cover the wavelength range 
5000 < A < 9200 A, galaxies at z = 3.5 contribute to 
the rest frame 1100 — 2050 A, while those a.t z — 4.5 
contribute to 900 — 1650 A. At longer rest-frame wave- 
lengths the composite spectrum is therefore largely con- 
tributed by galaxies at lower redshift. In addition there 
are wavelength-dependent sources of noise discussed in 
i j3.2l In summary, the data used to construct the com- 
posite spectrum in Figure [3] correspond to a mean red- 
shift z = 3.98 at = 1100 A, z = 3.95 at 1200 A, 
z = 3.90 at 1500 A, and z ^ 3.82 at 1650 A. For the 
wavelength range of interest in this work, the redshift 
bias Az < 0.15 is not particularly troublesome. 

3. COMPOSITE SPECTRA 

We now turn to the presentation of the composite spec- 
tra. We need to account for the difference between the 
redshifts determined using Lya on the one hand and the 
low-ionization interstellar lines on the other hand and 
the systemic redshift prior to co-addition. We also seek 
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FIG. 2: Redshifts and UV luminosities for the sample of 
131 galaxies with suitably accurate redshifts measured either 
from Lyaor strong interstellar absorption lines. Points are 
color-coded according to the equivalent width of Lya , WLya ■ 
Galaxies brighter than an apparent magnitude z'^^g < 26.0 
(solid line) typically have well-detected continua suitable for 
forming a composite spectrum. Galaxies with redshift z < 4.5 
(dashed line) are used for the composite spectra discussed in 
<J31-3ni while iJ5] includes galaxies at all redshifts shown here. 

to understand the signal/noise of the composite in terms 
of the statistical uncertainties and the variance among 
the individual spectra used to construct the composite. 

Composite spectra are constructed by shifting the indi- 
vidual spectra into the rest frame according to a deduced 
systemic redshift and then averaging the set. In general 
terms we will first identify a subsample based on their 
observable properties. Each spectrum in the sample is 
shifted to the adopted rest frame and interpolated to a 
common wavelength scale with a dispersion of 0.12 A. 
All spectra are normalized to have a median = 1 in 
the range 1250-1500 A. Spectra taken with DEIMOS are 
smoothed to a resolution of 1.9 Ato match the lower res- 
olution of F0RS2 data. The spectra are then averaged 
at each wavelength using a cr-clipped mean to reject out- 
liers arising from sky subtraction residuals and cosmetic 
defects. An equal number of positive and negative out- 
liers are rejected at each wavelength, totaling at most 
30% of the data. The remaining data are averaged with 
an arithmetic mean. 

Uncertainty in a composite spectrum will arise from 
both the finite signal to noise and the variance of 
the individual galaxies. For example, the variance in 
Lya equivalent widths in our sample is much greater than 
the uncertainty measured in the individual spectra. It is 
especially important to quantify the sample variance for 
weak features that are generally not detected in indi- 
vidual spectra. We account for sample variance with a 
bootstrap technique. For each composite spectrum we 
create 100 alternate composites using the same number 
of spectra but drawn at random from the parent sample. 
Each alternate has an average 63% of the sample repre- 
sented with 37% duplicates. Every measurement made 
on the composite spectrum is repeated for each of the 
100 alternates. We then take the measurement error to 
be the standard deviation of the 100 alternate measure- 
ments, which reflects both the sample variance and finite 



signal to noise. 

As discussed in §2.3 the most challenging issue is 
to determine the systemic redshift prior to shifting to 
the rest-frame . Her e we follow the approach used by 
IShaplev et al.l (|2003D . As a first approximation we use 
the value of ZLya (where available) to construct a com- 
posite spectrum. This enables us to locate the stel- 
lar photospheric line CiiiA1176 in the composite where 
we detect a velocity difference of —330 kms~^ with re- 
spect to Lya . We can thus infer that Lya emission in 
our sample is redshifted on average by AvLya = 330 
kms~^. In a similar fashion, stacking spectra using the 
redshift zjs based on low-ionization interstellar absorp- 
tion results in a detection of CiiiA1176 with a veloc- 
ity offset of -1-190 kms~^. For comparison, at z ~3 
Stcidel et al (2010) find {AvLya) = +445 kms^^ and 
{Avis) = —164 kms^^. To construct composite spectra, 
we use either zi^y^ shifted by —330 km s or z/5 shifted 
by +190 kms~^to approximate the systemic redshift of 
each galaxy. We use the Lya -based redshift when avail- 
able since it is typically determined with greater precision 
than Zjs- Figure [3] shows the composite spectrum of 81 
galaxies in our sample with 3.5 < z < 4.5 and apparent 
magnitude z'j^g < 26.0 using this method. 

3.1. Uncertainties in the Systemic Redshift 

A natural concern is the extent to which these applied 
shifts might vary within the sample used to make the 
composite. This can be estimated from observations of 
higher si gnal to noi s e from spectra taken at lower red- 
shift. .Steidel et al.l ()2010[ ) quantify the offset between 
zis, ZLya, and the systemic redshift ZHa in a sample 
of 89 galaxies at z ~ 2.3. They find velocity offsets 
Avis = -170 ± 115 kms~^ and AvLya = 485 ± 185 
kms~^relative to Ha. Assuming this is representative 
of our data, the uncertainty in our systemic redshift is 
therefore likely to be cr(f) ^ 150 kms^^. An upper 
limit on a(v) can be estimated from the width of spectral 
lines in the composite spectrum. In particular the rest- 
wavelength of the stellar line CiilA1176 in the composite 
provides a valuable measurement of the average offset 
from the systemic velocity and its width provides an up- 
per limit on the effective spectral resolution. We measure 
a systemic velocity of 21 ± 101 kms^^ in the composite 
spectrum (Figure[3]) and a FWHM = 520 kms~^ (decon- 
volved from the instrumental resolution ~ 450 kms^^). 
The uncertainty in the adopted redshifts about the true 
systemic stellar value is therefore < 520 kms~^ FWHM 
or equivalently (7{z) < 220 kms~^, comparable in fact to 
w hat was achiev e d for i ndividual spectra at lower redshift 
bv lSteidel et"all (I20T0I) . 

3.2. Error spectrum 

Figure [3] shows the composite spectrum of our sam- 
ple as well as the la error spectrum derived using the 
bootstrap technique discussed above. The error at each 
pixel is calculated as the standard deviation of all av- 
eraged data points (excluding outliers), divided by the 
square root of the number of data points. There are sev- 
eral wavelength-dependent factors contributing to the er- 
ror spectrum in addition to the finite signal-to-noise of 
individual spectra. One factor is the intrinsic sample 
variance, seen clearly as a noise spike at the position 
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FIG. 3: The composite spectrum of 81 galaxies in our sample with 3.5 < z < 4.5 and apparent magnitude z'j^g < 26. The 
effective mean redshift for the sample averaged over wavelength isJ = 3.9. The strongest spectral features are labelled. The gray 
filled region shows the ±lcr error at each point, determined from the scatter of individual spectra used to create the composite. 
The error spectrum peak at 1216 A is due to large scatter in the intrinsic distribution of Lya equivalent widths. The error is 
lowest at ~ 1300 — 1500 A where the continuum signal to noise ratio is ~ 30. The error increases at shorter wavelengths where 
the instrument throughput is lower, and at longer wavelengths where sky emission is much stronger. 



of Lya , and also evident for CliA1334 and other ab- 
sorption features. Another is the decreased instrument 
throughput at A < 6000 A, leading to higher noise at 
shorter wavelengths. Similarly, stronger sky line emis- 
sion causes increasing noise at longer wavelengths. Fi- 
nally, the number of contributing spectra peaks at rest- 
frame A ~ 1150 — 1350 A , with increased noise at higher 
and lower wavelengths where fewer spectra are included. 
Our observed wavelength range 5000 < A < 9200 A cor- 
responds to a rest-frame 1000 < A < 1800 A at the mean 
redshift z = 3.9 of our magnitude-limited sample. Ulti- 
mately we acheive a signal to noise ratio in the continuum 
of > 10 between the Lya line and rest-frame 1800 A and 
~ 5 in the Lya forest, with a peak S/N ~ 30 at 1350 A. 

4. FEATURES IN THE COMPOSITE SPECTRUM 

We now discuss the composite spectrum (Figure ^ in 
more detail, focusing on the strong spectral features at 
1215-1550 A where we have the best signal to noise. In 
this wavelength range we detect Lya, NvAA1239,1243, 
SiiiA1260, Siii*A1265, OiA1302-hSiiiA1304 (blended), 
Siii*A1309, C11AI334, SiivAA1394,1403, Sin A1527, 



Siii*A1533, and CivAA1548,1550 at high significance. 
Our z = 3.9 composite spectrum is very similar to the 
composite of z = 3 LBGs presented in iShaplev et al.l 
(200(|), which we show in Figure 2] for comparison. The 
absolute magnitude distribution of our sample is broadly 
similar to that at z = 3; both cover the range —22 < 
AIuY < —20. Our disc ussion below f o llows closely that 
presented originally bv IShaplev et al.l (|2003| 1 but we are 
also interested in whether there are differences seen over 
the redshift range 3 < z < 4.5. We will discuss these 
possible evolutionary trends in §5. 

4.1. Lyman Break Galaxies: a physical picture 

It is helpful to begin by describing a possible physical 
picture of LBGs based on ma ny analyses o f the exten- 
sive observations at z ~ 3 (see lShaplevll2011l for a recent 
review). Typical C* LBGs at z = 3 have ultraviolet half- 
light radii ~ 2.0 kpc, stellar masse s ^ 3 x 10^° Mf?), and 
star formation rates ^ 50 M(7)yr~^(Bo"uwens et al ll2004l : 
iFerguson et"al|[200l IShaplev et al.ll200lD . Theltar for- 
mation surface density is sufficient to drive "superwinds" 
of outflowing gas, similar to those seen in local galax- 
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FIG. 4: Composite spectrum of all 81 galaxies in our sample with z < 4.5 and apparent magnitude z' < 26, compared to the 
composite spectrum of 811 LBGs at z = 3 presented in lShaplev et alJ (|2003l ). The higher redshift sample has a much stronger 
Lya forest break, slightly redder UV spectral slope, slightly stronger Lya emission, and weaker absorption lines. The inset shows 
a zoom-in of the region from 1200 — 1600 A, which contains most of the absorption lines of interest in this paper. 

ies where T,sfb. > 0.1 Moyr^ikpc-^ ()Heckmaij [20021 ) . 
Indeed, blueshifted interstellar absorption lines confirm 



there are outflows with typical velocities of ~ 150 km s 
and in some cases as hig;h as 800 kms~^(|Shaplev et al.l 
[200llPettini et al.ll2002tlQuider et al. 2009*. '20101). How- 
ever, the physical origin o f these outflows remains unclear 
(e.g. iMurrav et al.|[20Tol) . It is thought these outflows 
produce an extended circumgalactic medium (CGM) of 
ejected material. Outflowing gas is found in both low and 
high ionization states (e.g. Sin and Siiv). Low-ionization 
transitions such as Sin and Cn are mostly associated 
with neutral hydrogen whereas high ionization lines oc- 
cur in the fully ionized component . Based on trends in 
the strength of the various species, IShaplev et all (l200l 
suggest a geometry in which discrete clou ds of neutral 
gas ar e embedded in a halo of ionized gas. iSteidel et al.l 
(j^OlO") show that the neutral and ionized CGM compo- 
nents both extend to radii of at least 125 kpc. 

4.2. Lya 

Lyais the most prominent and diverse feature in our 
individual spectra. The line originates from hydrogen 
recombination in H II regions photoionized by massive 
stars and, as these stars dominate the adjacent contin- 



uum, its intrinsic equivalent width should be within the 
range W^r,„ry = 100 — 200 A for nearly all stellar popu- 
lations (e.g. iForero- Romero et al.ll201ll ). Our composite 
spectrum reveals a complex line profile with strong ab- 
sorption extending blueward to v = —4000 kms^^, and 
redshifted emission to i; = +1000 km s~^ (Figure [5]) with 
a peak offset oi v = +330 km s^^ relative to the systemic 
velocity. The net equivalent width WLya = 21 ± 3 A is 
only < 20% of the expected intrinsic value. 

The form of this Lya profile has been readily under- 
stood in terms of the physical picture discussed in 
Lya emission produced at the systemic velocity can es- 
cape only along a line of sight free of neutral hydrogen 
or if it is shifted in velocity far from resonance. As pho- 
tons escape they encounter blueshifted clouds of partially 
neutral gas, which absorb and re-emit isotropically. Pho- 
tons backscattered from neutral clouds at small radii will 
appear redshifted and have a higher chance of escaping. 
Additionally, photons scattered from neutral gas at the 
edge of the CGM can escape into the ionized intergalac- 
tic medium and wil l be observed as blueshifted emission 
(jSteidel et al.ll2010l ). Indeed, Figure [S] shows such a weak 
blueshifted emission peak in the composite spectrum. On 
the other hand, the reduced equivalent width compared 
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to that expected intrinsically could be due to many ef- 
fects including dust extinction, scattering at large radii 
where emission falls outside the spectroscopic slits and a 
non-zero escape fraction fesc of ionizing radiation. Since 
these effects cannot easily be disentangled, the absolute 
strength o f the Lyg line ruust be interpreted with caution 
(Paper Il. lSchenker eFalllmih . 

We can estimate the spatial extent of the neu- 
tral CGM (i.e. the radius at which escaping Lyais 
last scattered) based on the strength of the diffuse 
blueshifted Lya emission. The spatial profile of extended 
Lyg emission has b een well-quantified at z ~ 2.6 by 
iSteidel et al.l ()201lD who find that diffuse haloes are a 
generic property of star-forming galaxies leading to a to- 
tal Lya flux ~ 5 x greater than that measured within the 
spectroscopic slits. 

Average Lya surface brightness profiles are well-fit at 
large radii by an exponential form, viz. 



^Lya (b) = T,oexp{b/bi) 



(1) 



where bi ~ 25 kpc. The blueshifted emission shown in 
Figure [5] has an equivalent width W = 1.5±0.3 A mea- 
sured within a slit aperture of ~ I'.'O x I'.'O, corresponding 
to 7.1 X 7.1 kpc^ at 2 = 4. Assuming this approximates 
the peak value, we estimate So = 1-5 ± 0.3 A per 50 
kpc^ and integrating Equation [1] yields a total equivalent 
width of Wfiaio = iirbfJ^Q. Using the measured value of 
Eq, the scale length is given by 



k = (2.3±0.2)JW^haio/Akpc 



(2) 



Although the total Lya flux is not measured, we can 
estimate its value from theoretical expectations as well as 
observations at lower redshift. The equivalent width in 
Figure [5] is = 21.0 A (excluding the blueshifted emis- 
sion component), thus we can write Whaio — Wtot — 21 
A. As discussed above, the observed Wtot is diminished 
by dust and the escape of ionizing radiation. Indeed, 
ISteidel et all (1201 ID measure Wtot = 17 - 93 A in vari- 
ous subsamples of their data, all lower than the expected 
intrinsic value, Wtot=100-200 A . Considering the com- 
posite DEIMOS spectra s hown in Figure [51 and assum ing 
an intrinsic Wtot = 135 A ()Forero-Romero et al.ll20lTI ) we 
derive upper limits of 6; < 20 ± 2 kpc for galaxies with 
WLya > 0, and bi < 17 ±2 kpc for those with WLya > 20 
A. This latter constraint is the mos t stringent, and also 
likely closest to the true value of bi. ISteidel et al.l (|2011D 
flnd that their sample of Lya emitters (LAEs, defined 
as WLya > 20 A) has the highest Wtot (= 93 A) and 
largest scale length (6; = 28.4 kpc) of any sub-sample 
that they analyze. Our constraint of bi < 17 kpc for the 
LAEs therefore suggests that the characteristic size of 
Lya haloes may be smaller at z = 4 than at z = 2 — 3. 
Due to the large inherent uncertainties, direct measure- 
ments of low surface brightness Lya emission are needed 
to confirm this possibility. 

4.3. Low-ionization metal transitions 

According to the physical picture in H4.1[ absorption 
in low-ionization transitions occurs in both the inter- 
stellar medium and outflowing clouds of cool gas. The 
former is at the systemic velocity while the latter is 
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FIG. 5: Top: composite spectrum normalized to continuum 
flux levels, showing the velocity structure of Lya . The line 
profile consists of a broad blueshifted absorption trough at 

V > —4000 kms~^and strong redshifted emission extending to 

V — -f 1000 kms~^. Bottom: velocity profile of Lya showing 
a blueshifted secondary peak. The composite spectra in this 
plot are constructed only from DEIMOS data, with spectral 
resolution R ~ 2000. The composites show a blueshifted 
Lya emission peak centered at ii = —600 kms~^, correspond- 
ing to the maximum outfiow velocity of neutral gas as seen in 
the average velocity profile of low-ionization absorption lines 
(black). Blueshifted Lya emission arises from photons scat- 
tered at the leading edge of outfiowing neutral CGM. 

blueshifted. Such transitions are gen erally saturated 
(jShaplev et al.ll2003l:lPettini et al.l20"0^ so the line depth 
at a given velocity provides a measure of the areal cov- 
ering fraction fc of O and B stars by neutral gas along 
the line of sight. Specifically, the line profile is given 
by f{v) = /o • (1 — fc(y)) where /o is the continuum 
flux. While typical LBGs are too faint for detailed line 
profiles, high-resolutio n spectra of a few bright lensed 
sources at z = 2 — 3 (jPettini et al.l 12002) : IQuider et al.l 
120091 120101 ) have revealed absorption velocities ranging 
from ^ —1000 to -1-500 kms""'^ with the highest cov- 
ering fraction at w ^ —200 kms~^. The mean low- 
ionization absorption line velocity offset in our composite 
is vlis — —190 kms~^in good agreement with t hat at 
z = 2-3 (|Shaplev et al.ll2003HSteidel et aLllMoh . 
Of particular diagnostic value are the two relatively 
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FIG. 6: (Left:) comparison of the absorption hne profiles for 
Sin transitions at 1260 and 1527 A in our composite spectrum. 
(Right:) The profile ratio (1 — /i26o)/(l — /1527) with the 
optically thick and thin regimes indicated; the shading refers 
to the la uncertainty. Although both transitions are optically 
thick at low velocities \v\ < 200 km s ~ ^ , there is some evidence 
of optically thin gas at \v\ > 200 kms~^. 

unblended transitions of Sin at 1260 and 1527 A whose 
equivalent width ratio is a valu able tracer of the opti- 
cal depth (jShaplev et all l2003i) . The profiles of both 
lines and their ratio are shown in Figure |6l Gas near 
the systemic velocity {\v\ < 200 kms""'^) is clearly opti- 
cally thick, while for v < —200 kms^^ it is intermediate 
in optical depth. This difference is also seen in high- 
resolution spectra of lensed galaxies ([Pettini et al.ll2002l : 
iQuider et al.ll2010[ ) , and optically thin g as is seen at large 
galactocentric radii (jSteidel et al.ll2010l ) suggesting that 
smaller, optically thin clouds are more easily accelerated 
to high velocity and large distances. We note that the 
column density at which both Sill transitions become 
optically thick (r > 1) is iVsin = 1-9 x 10^^ cm"^. As- 
suming log(Si ii/H) = —4.86 as m easured for the lensed 
galaxy cB58 (|Pettini et al.l 120021 ) . we estimate that op- 
tically thin absorption occurs in clouds with hydrogen 
column densities A'hi < 10^^ cm~^. 

4.3.1. Fine structure transitions 

A satisfying aspect of our composite spectrum is the 
successful identification of fine structure emission lines 
of Sin (see marked features in Figure [3]) . Sin ions in 
the CGM absorb photons in the resonance transitions 
and immediately re-emit a photon at approximately the 
same velocity. A photon absorbed at 1260 A will be 
re-emitted at either the same wavelength or to the fine 
structure transition Sin*A1265. Likewise a photon ab- 
sorbed at 1527 A will be re-emitted as either SinA1527 or 
Sin*A1533, and an absorbed SinA1304 photon can be re- 
emitted as Sin*A1309. In all three cases the probability 
of emission in the resonant and fine structure transitions 
is approximately equal. 

Since no absorption is seen in the fine structure tran- 
sitions, we infer that atoms in the excited ground state 
will typically decay to the ground state before absorb- 
ing another photon. Since every absorbed photon is re- 
emitted, the net equivalent width of the resonant and 



fine structure transitions is W^sin + Wgiu* = 0. The pre- 
cise equivalent widths depend on the initial absorption 
W^SiiLabs and the optical depth. If the gas is optically 
thin, re-emitted photons will immediately escape giv- 
ing W^Sin,em = W^Sin*,cm = -0.5W^siii,abs- In the limit 
of optically thick gas, resonant photons will be con- 
tinuously scattered until they emerge as Sin* (after 2 
scatterings on average). In this case, VFsin.om — and 
Wsiii*,om = ^W^Siii.abs- Sincc the majority of absorption 
occurs in optically thick gas ( ij4.3p . we expect the equiv- 
alent width of Sin absorption lines to reflect the kine- 
matics and covering fraction of neutral gas with minimal 
contamination from Sin re-emission. 

4.3.2. The spatial extent of low-iomzation absorption 

From the picture above ( §4.3.1|) . we expect that the 
equivalent width of the fine structure emission lines 
should be equal and opposite to the resonant line equiva- 
lent width. This is not the case: W^i265/W^i260 = — 0.56± 
0.36 and VK1533/VF1527 = -0.52 ± 0.20. (VF1309/VF1304 is 
more difficult to quantify since SinA1304 is blended with 
O1AI3O2). There are two possible explanations. One is 
that scattered photons have larger path lengths and are 
subject to greater dust attenuation. Since the UV con- 
tinuum slopes imply little differential extinction (mean 
E(B-V) — 0.10) and we expect only 2 scatterings for 
the average Sin* photon, this seems unlikely. More rea- 
sonably, the emitting region could be larger than that 
sampled by our slits. Although the I'.'O x I'.'O slit aper- 
ture samples most of the continuum light and its line of 
sight absorption, the CGM which absorbs and isotropi- 
cally r e-emits scattered photons extends to much larger 
radii (jSteidel et al.l [20la [20T1I ). The strength of fine 
structure emission provides a direct measurement of the 
amount of absorption at small radii contained within the 
slit. The fine structure to resonant absorption line ra- 
tio suggests that a fraction 0.53 ± 0.17 (combining both 
A1265 and A1527 measures) of the total fine structure 
emission is contained within the extraction aperture of 1 
arcsec^. The half-light radius of Sin* emission, and hence 
Sin absorption, thus corresponds to ^ 0'.'5 — 3.5 kpc at 
z = 4. This scale is only slightly larger than the median 
half-light radius rh = 0'.'31 of galaxies in our sample, and 
would be reached in only 20 Myr at the typical outflow 
velocity (190 kms~^). 

It is inst ructive to reconsider t he z = 3 composite spec- 
trum from'Shaplev et al.l (|2003l ) where the fine structure 
lines are clearly seen (but were not interpreted fully along 
the discussion above at the time). The Sin* lines are 
noticeably stronger in our z ~ 4 composite, with an av- 
erage Wsiu* = 0.7 ± 0.2 at z ~ 4 c.f. 0.30 ± 0.05 A at 
z = 3. Furthermore, the Sin absorption lines are weaker 
in the z = 4 composite, suggesting the absorption takes 
place at larger radii at z = 3. This difference cannot be 
explained through instrumental differences between the 
two surveys. Smoothing our z = 4 com p osite to match 
the 3.25 A resolution of iShaplev et al.l (|2003[ ) reduces 
the equivalent widths by < 5%. Likewise the wider slits 
used for the z = 3 data (I'.'4, 11 kpc) c.f. the z — 4 
data (I'.'O; 7 kpc) is not the cause. If the CGM prop- 
erties are similar, we expect the fine structure strength 
to constitute a larger fraction of the resonant absorption 
line strength in the z = 3 data, contrary to observations. 
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Defining Rps = -(W^i265 + VF1533)/ (M^i260 + ^^1527), we 
find Rfs = 0.53 ±0.17 for the z = 4 composite spectrum 
and Rfs = 0.16 ± 0.04 for z = 3, indicating a smaller 
characteristic radius of fine structure emission at z = 4. 

Inescapably, therefore, we conclude the circumgalac- 
tic gas around LBGs at z = 4 differs in two important 
ways from that at z = 3: there is greater low-ionization 
absorption at small radii at z = 4 and less total low- 
ionization absorption. In addressing the origin of this 
effect, Figure [TOl shows that the difference in equivalent 
width is at least partially linked to the kinematic offset 
from the Lya emission. Thus it seems a higher covering 
fraction at small radii is required to produce the stronger 
fine structure emission. Higher resolution observations 
of individual bright, or perhaps gravitationally-lensed, 
galaxies at z ~4 will ultimately be required to separate 
the relative contributions of covering fraction and kine- 
matics in explaining this result. 

4.4. High-ionization lines 

The high-ionization lines Siiv and Civ arise both in 
interstellar gas and in stellar P-Cygni winds. The veloc- 
ity centroid of Siiv, v — 140 kms~^, is consistent with 
the low-ionization interstellar absorption lines suggesting 
that most absorption is interstellar in origin. This is sup- 
ported by the absence of a significant redshifted emission 
component expected for a P-Cygni profile. In contrast. 
Civ is broader with a larger absorption velocity offset 
V = —370 kms~^and redshifted P-Cygni emission indi- 
cating a large contribution from stellar winds. We also 
detect the NvAA1239,1243 P-Cygni feature, ahhough the 
proximity to Lya makes this feature difficult to study in 
detail. 

We can determine the optical depth of highly ion- 
ized outflowing gas from the ratio of Siiv absorption 
lines. The ratio W^siivAi394/W^SiivAi403 = 2.0 for opti- 
cally thin abso r ption , and 1.0 in the optically thick case. 
IShaplev et al.l ()2003D find optically thin absorption in 
composite spectra of LBGs at z = 3, whereas we mea- 
sure a ratio 1.4 ±0.4 in the composite spectrum shown in 
Figure [3] indicating a significant contribution of optically 
thick absorption at z ~ 4. The total equivalent width of 
the Siiv doublet is weaker by a factor 0.72 ± 0.12 in Fig- 
ure [3] compared to the z = 3 composite of IShaplev et al.l 
([2003). The combination of higher optical depth and 
lower equivalent width suggests that the velocity range 
and/or covering fraction of the ionized gas traced by Siiv 
is lower at higher redshift. 

4.4.1. Metallicity 

The P-Cygni profile of CiV is sensitive to metallicity, 
and the combination of CiV and Hen equivalent widths 
constrains both the age and metallicity. Here we compare 
the equivalent widths measured in the composite spec- 
trum (Figure ^ with theoretical models in order to esti- 
mate the typical metallicity of galaxies in our sample. We 
note that the equivalent width of CiV contains significant 
interstellar absorption, so the value reported in Table [T] 
should be treated as an upper limit on the P-Cygni com- 
ponent. Assuming the interstellar absorption component 
of Civ is similar to that of Siiv {W ~ 1.0 A), we take the 
P-Cygni absorption component to have equivalent width 
— 1.6 ± 1.0 A with a conservative uncertainty. We com- 
pare this estimate and the measured equivalent width 



of Hell (reported in Table [T]) with predictions from the 
stellar population synthesi s code BPASS presented in 
lEldridge fc Stanwavi (|2009t ). We consider BPASS mod- 
els which include binary evolution with continuous star 
formation rate, and determine the difference AW be- 
tween observed and predicted equivalent width as a func- 
tion of metallicity and relative carbon abundance. We 
restrict the stellar population age to be within the la 
scatter of the median value for galaxies in our sample, de- 
termined to be lO^ '^iO-S years from spectral energy den- 
sity fit s assuming constant star formation with a Kroup^ 
(pO^) initial mass function. Figure [7] shows the result- 
ing relative error between BPASS models and measured 
equivalent width, which we define as 



Relative error = 0.5 



AW, 



Civ 



CTCiv 



/AT^Hcii 



CTHoii 



(3) 



Relative error values < 1 are thus consistent with the 
data. 

Models with Z < 0.004 and somewhat depleted carbon 
abundance are in good agreement with the data. We note 
that observations of both local and high-redshift galax- 
ies indicate typical relative carbon abundances (C/0) 
~ 2 — 5x lower than the solar value (equivalent to 
Xc = 0.2-0.5 in Figure [T]) for metallicities Z < 
0.004 (iKolD ulnickv & S kiThmaiT9 98: Shap lev et al1l2Q0l 
lErb et al.ll2010. : .Eldridge fc Stanwavi.201l[ ). Solar metal- 
licity models {Z = 0.020) do not fit the data. We con- 
clude that the typical metallicity of galaxies in our sam- 
ple is Z < 0.004 or < 0.2 x solar metaUicity. Bright 
galaxies at z < 3.8 in the F0RS2 sample have gas- 
phase metallicity 12 ± \og{0/H) = 7.7 — 8.5 or about 
0.1 — 0.6 X the sola r value, in reasonable agreement 
(jMaiolino et al.|[2008fl . These values are consistent with 
the allo wed metallicity range 0.1 < Z/Zq < 0.6 inferred 
for the IShaplev et a l. ( 2003 ) composite us ing the same 
BPASS models (lEldridge fc Stanwav l l20ll. Given that 
the eq uivalent widths of CiV and Hell in IShaplev et al.l 
(|200a agree well with those in Table [1] we expect the 
typical metallicity of galaxies to be similar in both sam- 
ples. 

5. SPECTROSCOPIC TRENDS 

We now turn to an analysis of how the various spectro- 
scopic features discussed in §4 are related to observable 
properties of LBGs as a prelude to considering how they 
might evolve with redshift. We will begin with the de- 
pendence of low-ionization absorption line strength with 
Lya equivalent width. This is motivated by the common 
physical dependence of these features - both are gov- 
erned by the kinematics and covering fraction of neutral 
circumgalactic gas - and also beca use this has be en ex- 
amined in detail for LBGs at z = 3 (jShaplev et al.l l2003'). 

We construct composite spectra of three sub-samples of 
galaxies divided according to their Lya equivalent width. 
Defining Wlis sls the average equivalent width of A1260, 
A1303, A1334, and A1527 A features. Figure M shows 
weaker Wlis for galaxies with s tronger Lya . This wa s 
also noted for z = 4 LBGs by iVanzella etldl ()2009D . 
Both the z=3 and z=4 samples show this trend, al- 
though Wlis is weaker at z = 4 for galaxies with 
strong Lya emission. This difference can be attributed 
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FIG. 7: Relative error in the equivalent widths of Civ 
and Hell predicted by the stellar population synthesis code 
BPASS and values measured for the composite spectrum 
shown in Figure [3] (see text for details). Relative error val- 
ues < 1 indicate good agreement. Model results for a range 
of metallicity Z and carbo n depletion factor X c as a func- 
tion of age are described in lEldridee &: Stanwavl (j2011h . The 
values Xc = 0.5 (0.1) correspond to a carbon abundance re- 
duced by a factor of 2 (10) relative to solar abundance ratios. 
Lines are plotted for the median age of our sample (^iQ^-^^os 
years) determined from spectral energy density models with 
the same initial mass function and star formation history used 
in BPASS, and error bars indicate the allowed range for ages 
within la of the median. Measurements of C/O abundance at 
both low and high redshift indicate carbon depleti on factors 
Xc ^ 0.2 - 0.5 in galaxies with Z < 0.2 Zq (e.g. lErb et all 
I2010f ). For this range of Xc, models with Z < 0.004 (equiv- 
alent to Z = 0.2 Zq) are in good agreement with the data. 
Solar metallicity models are unable to reproduce the observed 
equivalent widths. 

to the luminosity-dependent trend of stron ger Wlvoi and 
weaker Wtjfi in fai nter galaxies (Paper I; iShaplev et al.l 
[T003; Van zella et a l. 2009). If we consider a subset of 
galaxies in our sample with similar absolute magnitudes 
to those observed at z = 3 (—21.5 < Muv < —21.0), we 
recover the same normalization (Figure [SJ . Our sample 
is 90% complete at the corresponding apparent magni- 
tudes, so we expect a negligible bias. 

Having established the correlation of Wlis with 
WLya , we now examine the dependence with other de- 
mographic properties. We divide the full sample into two 
bins of equal size according to each property of interest. 
The results are shown in Figure [5] We briefly review the 
trend of each property with Wlis and WLya and dis- 
cu ss the physical origin . Man y of these trends were seen 
in IShaplev et all ([2003!) and IVanzella et all ([2009) and 
are clearly inter-related due to correlations between the 
demographic properties. 

Less luminous galaxies have stronger W^ya and weaker 
Wlis illustrating that higher star formation rates drive 
larger amounts of neutral gas into the CGM with higher 
velocity and/or covering fraction. Defining the ultravi- 
olet spec tral slope as = 5.30 • (z^b — zab) — 2.04 for 
B-drops (|Bouwens et al.|[2009[ ). we also find that bluer 
galaxies with lower /3 have stronger WLya and weaker 
Wlis- This trend was also noted in Paper I, which 
showed that LBGs with strong Lya emission have sys- 



tematically bluer /3. Since neutral gas also presumably 
contains dust, the same gas which gives rise to Wlis 
also reddens the continuum. Stellar masses are mea- 
sured for 60% of our sample for which there is unconfused 
Spitzer/IRAC photometry. As expected from the trends 
with luminosity, lower mass galaxies have stronger WLya 
and weaker Wlis- Finally, measuring half-light radii r/,, 
from the GOODS ACS data with Sextractor, we find 
smaller galaxies have stronger WLya and weaker Wlis- 

Could these trends could be due to selection effects or 
sample bias? The trend of stronger Lya for less luminous 
LBGs is of particular concern, since fainter galaxies will 
require a larger WLya for detection. To address this, we 
consider only the 32 galaxies with apparent magnitudes 
z'j^g < 25 for which the spectroscopic sample is 90% com- 
plete. This sub-sample is divided into two equal bins and 
the results confirm that trends seen in the larger sample 
also hold in brighter galaxies unaffected by sample bias. 

The trends shown in Figure [S] are generally consistent 
with the overall trend of Wlis with WLya (Figure [H) . 
The WLya variations mostly arise from the distribution 
of neutral gas (traced by Wlis), with relatively little ef- 
fect from other demographic properties examined {Muv, 
f3, M*, r/i). However, demographic properties do have 
some effect. The strongest deviation seen in Figure [S] 
is that with (3, which shows stronger Wlis than would 
be expected at a giv en W^va fo r galaxies with red UV 
slopes. Noted also bv IShaplev et al. (2003), this suggests 
that outflowing neutral gas contains dust which reddens 
the continuum. Also, more luminous galaxies (i.e. those 
with higher star formation rates) have stronger Wlis 
at a given WLya (Figure E]). This is likely due to in- 
creased absorption at large velocit ies in galaxies wit h 
higher SFR, as observed at z ~ 1.4 (jWeiner et al.ll2009( ). 
For all other demographics, the composite spectra have 
values of Wlis within Icr of that expected purely based 

on WLya - 

In summary, the trends seen in Figure [5] arise almost 
entirely because of variations in the neutral gas covering 
fraction and/or kinematics, which are themselves corre- 
lated with the demographic properties. 

5.1. Kinematics 

We have now established that outflowing neutral gas 
is the dominant factor in determining both WLya and 
Wlis- However, variations from the trend of Wlis 
with WLya are apparent, particularly with (3 and Muv 
as discussed above. The most obvious mechanism for 
this behaviour is a systematic difference in the cover- 
ing fraction fc and kinematics of neutral gas. Assuming 
Wlis J fc{v)dv (where v is the outflow velocity), a 
higher fc and lower velocity range can conspire to give 
a constant Wlis- However, changing fc and v will also 
affect the transmission of Lya photons resulting in a dif- 
ferent WLya - It is therefore of interest to consider how to 
distinguish between the covering fraction and kinematics 
of neutral gas. 

The kinematics of neutral gas can be roughly param- 
eterized by the velocity dispersion and centroid of low- 
ionization interstellar absorption lines. We have mea- 
sure the velocity Av of low-ionization absorption lines 
with respect to Lya as a proxy for outflow velocity. Av 
is strongly correlated with WLya and Wlis in LBGs at 
z = 3 ((Shaplev et al..,2003i) in the sense that larger ve- 
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• z = 3 (Shapley et ol 2003) 

♦ z = 4 (this work) 
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FIG. 8: Equivalent width of low-ionization absorption lines 
compared to that of Lya . Gray diamonds are from our sam- 
ple at 2 ~ 4 divided into bins of WLya < 0, WLya =0 — 40 
A, and Whya > 40 A. Black circles show the sample of Shap- 
ley et al (2003) at mean z = 3, divided into quartiles of 
WLya ■ Galaxies in the z — 3 sample have typical luminosities 
—21.5 < Muv < —21.0 corresponding to the blue triangles. 
Galaxies of the same luminosity lie on the same correlation 
between WLya and Wlis at both z = 3.9 and 2 = 3. Fainter 
galaxies have weaker low-ionization absorption lines at fixed 

WLya ■ 
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FIG. 9: Equivalent width of low-ionization absorption lines 
compared to that of Lya , divided according to observable de- 
mographic properties as described in the text. Points binned 
by WLya at z = 3 and 2 = 4 are the same as in Figure [S] The 
sample was divided into two bins for each demographic prop- 
erty {Muv, M,, P, and r^), with WLya and Wlis measured 
from composite spectra of the galaxies in each bin. 

locities are associated with stronger interstellar absorp- 
tion and weaker Lya emission. Figure [TU] shows Av mea- 
sured from the same composite spectra used to deter- 
mine demog raphic trends ( e .g. Fi gure Et, as well as the 
resuhs from iShaplev et "all (|2003l ). Ah z ~ 4 compos- 
ites are consistent (within la) with the relation mea- 
sured at z = 3 as we ll as the mean Av m easured for 
B-dropout galaxies by iVanzella et all (|2009[ ). We mea- 
sure a velocity dispersion for each composite, and find 
that each is within 1.2ct of the effective spectral reso- 



FIG. 10: Kinematic offset Av between Lya emission and 
low-ionization absorption lines, as a function of WLya ■ The 
composites used to measure At; and Wlvo are the sam e as 
in Figure [51 Data at z — 3 are from Sha pley et al.l (|2003l ). 

lution a = 290 km s^^ (measured for the stellar [Cm] 
feature in Figure [3]). We are therefore unable to detect 
trends in outflow kinematics with demographic proper- 
ties in composite spectra. Higher signal-to-noise data, 
higher spectral resolution, or detailed studies of individ- 
ual galaxies are required to address trends in kinematics 
and covering fraction of neutral gas at z = 4. 

6. THE EVOLVING CGM 

We now turn to a discussion of the redshift evolution of 
the neutral CGM surrounding typical LBGs. The most 
useful probes are Lya and low-ionization absorption lines 
which trace the kinematics and covering fraction of neu- 
tral gas, and fine structure emission lines which provide 
a constraint on the spatial extent of the absorbing gas 

(SXH). 

The common dependence of Lya and low-ionization ab- 
sorption lines on the neutral CGM results in a strong 
correlation between Wlis smd WLya ■ Various physical 
properties of LBGs are correlated with both Wlis and 
WLya , but in such a way that the relation between Wlis 
and WLya remains nearly constant (Figure [H])- Further- 
more, the WLis-WLya relation at fixed Muv does not 
change significantly with redshift between z = 3 and 
z = 4 (Figured]). 

To further examine evolutionary trends with redshift, 
we now divide our spectroscopic sample into two bins 
of redshift at fixed Mjjv, now including galaxies at all 
redshifts (no longer restricted to z < 4.5 as in previ- 
ous sections). We consider galaxies with absolute mag- 
nitude —21.5 < Mrr v < —20.5, cho s en to be representa- 
tive of the sample in IShaplev et al.l (|2003f ). There are 64 
galaxies in our sample within this Muv range (see Fig- 
ure H]) . We construct composite spectra of galaxies with 
redshift above and below the median z — A.l and mea- 
sure the equivalent width of Lya and low-ionization lines 
(both resonant absorption and fine structure emission) in 
each composite. The results are given in Table [2] along 
with the demographic properties of galaxies in each sub- 
sample. The quantities Muv, and rh are measured 
from photometry while /S is determined from a direct fit 
to the ultraviolet continuum in the composite spectrum. 
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We define fx oc A'^ and fit the rest frame 1300 - 1700 
A to determine /3, with uncertainty quantified using the 
bootstrap method described in ij3l Aside from redshift, 
the demographics of each sub-sample are quite similar. 
The higher redshift galaxies have slightly higher average 
M,, smaller r^, and smaller (bluer) /3. WLya is consis- 
tent for both to within the sample variance, and is also 
consistent with the value WLya = 14.3 A measured for 
the composite spectrum of z = 3 LBGs in Shapley et al 
(2003). The most striking difference is in the strength 
of the low-ionization absorption lines, which are signifi- 
cantly weaker at higher redshifts (Figure [TT|) . The vari- 
ation in Wlis is ^oi explained by systematic differences 
in WLya or demographic properties, hence we seek an 
alternate explanation. 

We first examine whether the evolution in Wlis could 
arise as a result of different equivalent width distribu- 
tions for Lya . Although the mean WLya across our two 
redshift subsamples is similar, the lower redshift sub- 
sample has a broader distribution and contains more 
galaxies with Lya in absorption (WLya < 0). This is 
reflected in the larger sample variance in WLya at lower 
redshift (Table [2]) . We evaluate the effect of this poten- 
tial bias on Wlis by constructing a composite spectrum 
from a subset of the z < 4.1 galaxies with intermedi- 
ate WLya — — 30 A, resulting in a consistent mean 
WLya with sample variance reduced by a factor of 2.5. 
This composite has Wlis = —1-5 A, 0.2 A higher than 
when the full range of WLya is used, but still consider- 
ably lower than the value Wlis = —1-0 A measured for 
the higher redshift galaxies. We therefore conclude that 
differences in the WLya distribution are insufficient to ex- 
plain the observed variation in absorption line strength 
with redshift. 

There are several possible physical explanations for 
the evolution of Wlis with redshift shown in Figure fTTl 
One possibility is that the kinematics and/or covering 
fraction of neutral gas are systematically different. For 
example, an outflowing wind with fixed input energy 
and momentum will reach higher velocity at lower red- 
shifts due to lower density of the IGM. However, this 
effect should be stronger between z = 3.0 — 3.8 than 
from z = 3.8 — 4.7, whereas the decrement in Wlis is 
much stronger from z = 3.8 — 4.7 (Figure fTT|) . Further- 
more, we measure a higher offset between the velocity 
centroid of Lya emission and low-ionization absorption 
for z > 4.1 galaxies (Aw = 660 ± 80) than for z < 4.1 
{Av — 550 ± 40). Both measurements are consistent 
with the trend shown in Figure [10 Based on the tren d 
of Wlis with Aw seen at z = s lshaplev et al.ll2003f ). 
we would then expect stronger Wlis at higher redshift. 
Kinematics are thus unable to explain the difference in 
Wlis, at least with the information currently available. 
Data with higher spectral resolution and signal-to-noise 
are needed to fully address differences in the covering 
fraction and kinematics of neutral gas. 

Another obvious potential cause of weaker Wlis, or 
equivalently decreased WLya at a fixed Wlis, is the 
presence of neutral hydrogen with very low column den- 
sity of heavy elements. In this case we would expect 
optically thin absorption line profiles. This could be di- 
rectly tested as shown in §4.3.1) but the signal-to-noise 
of SillA1527 in the high redshift composite is too low 



to constrain the optical depth. We do see some evi- 
dence that the column density of neutral gas is lower 
at higher redshifts based on Lya line profiles in the com- 
posite spectra, shown in Figure 1121 The Lya absorption 
trough at - 1200 - 1215 A is significantly weaker in the 
higher redshift composite. This absorption arises at least 
in part from damping wings of high column density gas 
with A'hi > 10^" cm~^ associated with the low-ionization 
metal absorption lines (e.g. iPettini et aIll2000L [20021. 
The higher redshift galaxies are therefore characterized 
by lower typical iVui or/and a lower covering factor of 
high column density gas. Higher signal-to-noise data at 
z ~ 5 is required to determine whether this affects the 
optical depth of low-ionization metal transitions. 

Weaker low-ionization lines could also be caused by a 
systematically higher ionization state at higher redshifts 
resulting in lower column densities of low-ionization 
gas. This would also explain the weaker damped 
Lya absorption trough seen at higher redshifts (Fig- 
ure I12[) . This scenario can be tested by measuring the 
equivalent widths of higher-ionization silicon transitions, 
in particular SiiiiA1206 and SiivAA1393,1402. We find 
that the Siiv lines are weaker at higher redshift by a fac- 
tor of 0.68 ± 0.25 in VFsiivAAi393,i402, consistent with the 
decrement in Wsiu- We do not detect the SiiiiA1206 line 
in the high redshift composite, with a la upper limit of 
W^Sim < 0.8 times that of the lower redshift composite. 
SiiiiA1206 is marked in Figure [12] and is clearly stronger 
at lower redshift. We therefore find no evidence of a sig- 
nificant change in the ionization state of outfiowing gas 
at higher redshifts. 

A final possibility for lower column density of heavy 
elements is that Lya is scattered by "cold streams" 
of nearly metal-free neutral gas, predicted by simula- 
tions to accrete onto galaxies at small ra dii and with 
increasing rates at higher redshi ft (e.g. iDekel et"al] 
120091: iFaiacher-Giguere et al.|[201lD . If they exist, such 
streams could be identified as Hi absorption systems 
with no correspond ing absorption from heavy elements. 
iSteidel et all ()2010f ) find little evidence for the presence 
of such streams around LBGs at z = 2 — 3. A similar 
study at higher redshift would be very interesting but 
also extremely challenging with current observational fa- 
cilities. 

Could the difference in Wlis with redshift be caused 
in part by a different physical extent of circumgalactic 
gas? In ti4.3.2l we showed that fine structure emission 
lines can be used to measure the amount of absorption 
taking place within the size of the spectroscopic aper- 
ture. This quantity does appear to change with redshift, 
in the sense that absorption takes place at smaller char- 
acteristic radius in LBGs at higher redshift. The com- 
posite spectra of LBGs at z = 4 has weaker absorption 
lines and str onger fine struc ture emission than the z = 3 
composite of iShaplev et al.l {2003) (Figure llt, in qualita- 
tive a greement with trends seen at z = 3 (jShaplev et al.l 
l2003t see their Figure 9). We have only a weak con- 
straint on the redshift evolution at z > 4 due to limited 
signal-to-noise in the high redshift composite. The ratio 
of total Sin fine structure emission W^sin- in the high- 
and low-redshift composites is 1.1 ± 0.5 (taken as 



H^1265 + W'i309 + Wi533 
Wi2eo + Wi303 + W^1527 
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FIG. 11: Equivalent width of low-ionization absorption lines 
measured from composite spectra of LBGs at different red- 
shifts. The DEIMOS and FORS2 data presented in this paper 
are separated into two sub-samples of equal size as described 
in the text. We also show the equivalent resu lt at z = 3 
from the composite spectrum of Shaplc v ct al] (120031 ). with 
redshift distribution described in Stei del et al.l (|2003l ). The 
average low-ionization line profile of each composite is shown 
in the inset. Low-ionization absorption lines are significantly 
weaker for galaxies in the highest redshift composite. All 
three samples have consistent mean luminosity and WLya ■ 

with values given in Table[2]), while the corresponding ra- 
tio of absorption line strength Wlis is 0.6 ± 0.2. The in- 
ferred fraction of low-ionization absorption within the slit 
aperture is thus a factor 1.8± 1.0 larger in the higher red- 
shift sample, consistent with no evolution. Higher signal- 
to-noise is required to accurately constrain the spatial 
extent of low-ionization absorption at z > 4. 

To summarize, we find weaker low-ionization absorp- 
tion lines in LBGs at higher redshift with fixed Mjjv- 
The difference in Wlis is not consistent with the trends 
observed for Lya and various demographic properties ex- 
amined in |j5l We have discussed several possible causes 
for this discrepancy including variation in the kinemat- 
ics, covering fraction, optical depth, ionization state, and 
spatial extent of absorbing gas. The data show that the 
ionization state has no significant effect, but we are un- 
able to conclusively address other possible causes. Spec- 
tra of galaxies at z > 4 with higher signal-to-noise are 
required to more accurately constrain the optical depth 
and spatial extent. Additionally, spectra of individual 
galaxies taken with higher spectral resolution will be re- 
quired to independently determine the covering fraction 
and kinematic structure of absorbing gas. 

6.1. Galaxy evolution 

We argue in ^that Lya equivalent width is determined 
primarily by the neutral CGM, which is correlated with 
various demographic galaxy properties (Figure |9]) . It is 
well established that these demographic properties vary 
with redshift and we expect WLya to vary accordingly. 
Large photometric surveys have shown that Lyman break 
galaxies at increasi ngly higher re dshifts z > 3 have 
lower luminosities (iBouwens et al.l [2 007. 2011a|). bluer 
UV spectral slopes (Bouwens et al."2009, 20 llbl7 smaller 
stellar masses (,Stark et al..,2009; .Gonzalez et al...2011il . 



FIG. 12: Composite spectra of LBGs showing the region 
around Lya for the two redshift ranges described in SjB] We 
show the position of SiinA1206 blueshifted by 200 kms~^, 
Lya emission redshifted by 400 kms~^, and Lya absorption 
blueshifted by 1500 kms~^ where it is seen most prominently. 
The LyQ absorption trough is weaker in the higher redshift 
composite indicating a lower incidence of neutral gas with 
high column density. Sim absorption is also weaker in the 
higher redshift composite. 

and smaller sizes (jBouwens et al.l 120041 : iFerguson et al.l 
l2004f ). This is in accordance with inside-out galaxy 
growth: galaxies increase in size and stellar mass as they 
evolve with time, while increasing metallicity and dust 
content reddens the ultraviolet continuum. Simultane- 
ously, star formation drives large-scale outflows of gas 
which reach larger distances and are accelerated to larger 
velocities at later times (e.g. iMurrav et al.ll2010l ). Galax- 
ies which are more evolved (i.e. larger, more massive, 
redder) should therefore have a CGM characterized by 
larger spatial extent, larger velocity range, and higher 
covering fraction of neutral gas. Observationally this 
results in weaker WLya , stronger Wlis, and relatively 
weaker fine structure emission. These are pr ecisely the 
trends observed at both z = 3 ([Shaplev et al.|[2003.,') and 
z = 4 (this work). 

6.2. Lya in the epoch of reionization 

We reiterate that Lyman break galaxies at increasingly 
higher redshifts z > 3 have lower luminosities, bluer UV 
spectral slopes, smaller stellar masses, and smaller sizes. 
Notably, all trends in the demographics of galaxies at 
higher redshift are correlated with stronger WLya and 
weaker Wlis (Figure [Q]) . We therefore expect typical 
galaxies at higher redshifts to have, on average, stronger 
Lya emission. Earlier results from this survey confirm 
that strong Lya emission is more frequent in galaxies 
at higher redshift (Paper I; Paper II). We find no ev- 
idence of this trend reversing. In fact, galaxies with 
extremely small size, low mass, and blue /? tend to be 
the strongest Lya emitters (e.g. lErb et al] l2010t). We 
do, however, expect the average Lya emission strength to 
decrease significantly at increasing redshifts in the epoch 
of reionization due to ne utral hydrogen in the IGM (e.g. 
iHaiman fc Spaan£lll999f ). 

We have devoted considerable discussion to the prop- 
erties of Lya in part because Lya is of great interest as a 
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tracer of cosmic reionization. Several authors have now 
presented evidence that reionization was incomplete at 
z ~ 7 based on a rapidly decreasing fractio n of galax- 
ies with strong Lya emission at z > 6.5 (Sch enker et alJ 
[20Tlt lOno et aUlMlL IPentericci et alll201l[ ). Although 
trends at z = 3 — 4 suggest that the galaxies observed 
by these authors should have a higher fraction of strong 
Lya emission, we have presented evidence that the rela- 
tion between WLya and Wlis is systematically different 
at z > 4 (Figure [11]) . This implies a systematic difference 
in the spatial, kinematic, or optical depth structure of 
neutral circumgalactic gas compared to galaxies at lower 
redshift. The physical origin of this evolution and its 
effect on Whya at z > 4 will need to be understood 
in order to fully interpret the results of Lya surveys at 
higher redshifts in the context of reionization. 

7. SUMMARY 

The rest-frame ultraviolet spectra of star-forming 
galaxies contains a wealth of information about the prop- 
erties of the circumgalactic medium. In this paper we 
have presented an analysis of several features which trace 
the CGM with a focus on the properties of neutral gas. 
We find that the tre nds observed at lower redshift (z = 3; 
iShaplev et al.l [20031 ) also hold at z = 4 with approxi- 
mately the same normalization (Figure [8]). However, we 
find evidence for rapid evolution at z > 4 with lower 
Wlis at fixed Whya and Mjjv, suggesting a system- 
atic difference in the spatial distribution, kinematics, 
ionization state, or optical depth of circumgalactic gas 
at higher redshifts. We determine that the ionization 
state is not responsible for the observed evolution but are 
unable to distinguish between kinematics, covering frac- 
tion, optical depth, or the spatial extent of neutral gas as 
the likely cause. We are collecting additional spectra of 
LBGs at z > 4 with our ongoing survey, including high 
spectral resolution observations of bright lensed galaxies 
from which we can disentangle the kinematic profile and 
covering fraction of neutral gas. These data will allow us 
to address the precise magnitude and physical origin of 
evolution in circumgalactic gas properties. 

As a final note, we emphasize that the neutral CGM 
is of great interest in the context of reionization of the 
universe. Neutral gas in the circumgalactic medium ab- 



sorbs ionizing radiation, thereby inhibiting the ability 
of galaxies to reionize the universe. The escape frac- 
tion of ionizing photons is one of the most important 
and uncertain factors in determining the contribution 
of sta r-forming galaxies to reionization (Robertson et all 
l2010f l. We have shown that typical galaxies at higher red- 
shift have weaker low-ionization absorption lines based 
on their demographic trends, and presented new evi- 
dence that absorption lines are systematically weaker at 
z > 4 even for fixed demographic properties. This is 
likely caused by a lower covering fraction and/or veloc- 
ity range of neutral gas, and we will address the physical 
origin of this evolution with future data from our ongo- 
ing survey. Determining the redshift evolution of neutral 
gas covering fraction in LBGs will be of great interest 
for interpreting surveys of Lya emission in the context 
of reionization and addressing the role of star-forming 
galaxies in reionizing the universe. 
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Siiv 
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-0.8 ± 0.2 
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1526.71 
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1640.40 
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TABLE 1: Equivalent width and velocity of absorption and 
emission lines in the composite spectrum (Figure [3}. 



Property 


2 < 4.1 


2 > 4.1 


2 


3.76 ±0.21 


4.70 ±0.58 


Muv 


-21.0 ±0.3 


-21.0 ±0.3 


logM./M© 


9.6 ±0.6 


9.8 ±0.6 


rh (kpc) 


2.42 ±0.97 


1.94 ±0.74 




-2.02 ± 0.08 


-2.12 ±0.19 


WLyc. (A) 


15.4 ±5.9 


12.2 ±3.7 


Wl260 (A) 


-1.7 ±0.4 


-0.7 ±0.4 


1^1265 (A) 


0.8 ±0.3 


1.0 ± 1.0 


VVisos (A) 


-2.1 ±0.4 


-0.9 ±0.4 


Wl309 (A) 


1.0 ±0.3 


0.9 ±0.4 


Wl334 (A) 


-1.7 ±0.3 


-1.7±0.9 


Wl527 (A) 


-1.2 ± 0.3 


-0.6 ± 1.0 


VK1533 (A) 


0.7 ±0.3 


0.9 ±0.5 



TABLE 2: Mean demographic and spectroscopic proper- 
ties of LBGs at different redshifts. Error bars correspond 
to the standard deviation of values for individual galaxies in 
each sub-sample. /3 and equivalent widths of Lya and low- 
ionization metal transitions are measured directly from com- 
posite spectra, with error bars determined from a bootstrap 
method (see text for details). 



